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We have performed high precision measurements of the coherent neutron scattering lengths of 
gas phase molecular hydrogen and deuterium using neutron interferometry. After correcting for 
molecular binding and multiple scattering from the molecule, we find b np = (—3.7384 ± 0.0020) fm 
and b n d = (6.6649 ± 0.0040) fm. Our results are in agreement with the world average of previous 
measurements, b np = (-3.7410 ± 0.0010) fm and b nd = (6.6727 ± 0.0045) fm. The new world 
averages for the n-p and n-d coherent scattering lengths, including our new results, are b np = 
(-3.7405 ± 0.0009) fm and b ni = (6.6683 ± 0.0030) fm. We compare b nd with the calculations 
of the doublet and quartet scattering lengths of several nucleon-nucleon (NN) potential models 
and show that almost all known calculations are in disagreement with the precisely-measured linear 
combination corresponding to the coherent scattering length. Combining the world data on b n d with 
the modern high-precision theoretical calculations of the quartet n-d scattering lengths recently 
summarized by Friar et al, we deduce a new value for the doublet scattering length of 2 a n d = 
(0.645 ± 0.003(exp)±0.007(theory)) fm. This value is a factor of 4 more precise than the previously 
accepted value of 2 a„d = (0.65 ± 0.04(exp)) fm. The current state of knowledge of scattering 
lengths in the related p-d system, ideas for improving by a factor of five the accuracy of the b„ p 
and b n d measurements using neutron interferometry, and possibilities for further improvement of 
our knowledge of the coherent neutron scattering lengths of 3 H, 3 He, and 4 He are discussed. 

PACS numbers: 03.75.Dg, 07.60.Ly, 61.12.-q 

Keywords: neutron interferometry, scattering length, neutron optics, NN potentials, three- nucleon force, 
effective field theory, n-p, n-d, deuterium 



I. INTRODUCTION 

The three-nucleon system is both fascinating and re- 
markable. Despite decades of intensive study with in- 
creasingly sophisticated theoretical tools, fundamental 
facets of the system remain mysterious, confusing, and 
contradictory. Realistic nucleon-nucleon (NN) potentials 
underbind 3 He and 3 H by several hundred keV The p- 
d and n-d scattering lengths in the doublet s-wave chan- 
nel, which are apparently strongly correlated to the 3 He 
and 3 H binding energies respectively, are likewise poorly 
predicted by NN force models. Convergence between the- 
ory and experiment for these fundamental parameters 
can be obtained only by the ad hoc admixture of three- 
nucleon (3N) forces. Although it is well understood that 
3N forces must exist with a weaker strength and shorter 
range than the NN force, little else is known. The incom- 
plete nature of current 3N force models is demonstrated 
by their tendency to resolve certain problems, such as 
tri-nucleon underbinding, at the expense of exacerbating 
discrepancies in other observables 0,|ESIEIEIEIEI3- 

Of course one must first understand two-nucleon forces 
in detail before one can isolate possible 3N force effects. 
A wide variety of NN forces are deployed in modern 



nuclear physics calculations, but they have certain fea- 
tures in common. They all employ short-range, semi- 
phenomenological forces matched to a one pion exchange 
potential (OPEP) tail beyond w 1.4 fm, and they all give 
roughly the same results for the low-energy observables. 
They are all constrained to reproduce the NN observ- 
ables as well as the properties of the deuteron. They 
also must incorporate deviations from isospin symme- 
try to describe the data. One of the most sophisticated 
of the modern potentials-the AV18 potential-includes, 
in addition to purely electromagnetic terms, terms ac- 
counting for charge independence breaking (CIB) and 
charge symmetry breaking (CSB), which are phenomeno- 
logically adjusted to replicate the n-n, p-p and n-p scat- 
tering lengths [iH]. 

The accessible places to look for 3N force effects are the 
bound states of 3 H and 3 He and the scattering states of 
n-d and p-d. Since 3 H and n-d are free from electromag- 
netic complications (both theoretical and experimental), 
they are the systems of choice for precision tests. For 
all of these systems, the computational tools presently 
available are believed to be excellent. "Exact" Fadeev 
solutions are available for the scattering states below the 
deuteron breakup threshold. Above breakup threshold, 
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the n-d equations can also be solved exactly throughout 
the entire range of applicability of the potential mod- 
els employed. In addition to these exact models, other 
methods based on the reactance matrix (if -matrix), such 
as the hyperspherical harmonic expansion of Kicvsky et 
a2.[H|, the effective field theory of Bedaque et a/.^3], as 
well as dispersion theory/phenomenological models such 
as that of Hale|l3(, provide complementary insights and 
information that goes beyond tests of potential models. 

Recently, significant insight into certain features of 
three-nucleon systems has come from the effective field 
theory (EFT) approach 0G]1. EFT has been used to 
solve the three nucleon problem with short-range interac- 
tions in a systematic expansion in the small momentum 
region set by kb < I, where k is the momentum transfer 
and b is the scattering length. For the two-body system 
EFT is equivalent to effective range t heory and repro- 
duces its well-known results EUHlll. The chiral EFT 
expansion does not require the introduction of an oper- 
ator corresponding to a 3N force until next-to-next-to 
leading order (NNLO) in the expansion, and at this or- 
der it requires only two low energy constants 19]. With 
these two parameters determined from experiment, chiral 
EFT can make precise predictions for other three-body 
observables. It is clear that precision measurements of 
the three-nucleon scattering lengths provide both valu- 
able theoretical benchmarks and critical input parame- 
ters for the study of three-nucleon systems. 

The most important data, of course, is the measure- 
ment of the zero-energy scattering. With the neutron 
interferometer at NIST, we have made measurements of 
the coherent n-p and n-d bound scattering lengths, 6, 
which are related to the free scattering length a by 



Here, m is the mass of the neutron and M is the mass of 
the atom. For hydrogen, a is the linear combination of 
the singlet and the triplet scattering lengths given by, 

(1/4) l a np + (3/4) 3 a„ p , (2) 

and for deuterium it is the linear combination of the dou- 
blet and quartet scattering lengths, 

a nd = (1/3) 2 a nd + (2/3) 4 a nd . (3) 

These new measurements have uncertainties that are 
comparable to (in n-p) and smaller than (in n-d) pre- 
vious results. 

The n-d measurement yields a linear combination of 
the doublet and quartet scattering lengths and not the 
separate channel scattering lengths. There are good rea- 
sons however to have confidence that the doublet ( 2 S!/ 2 ) 
scattering length can be extracted from this measurement 
reliably with some help from theory. In the first place, 
the quartet s-wave scattering length ( 4 S 3 / 2 ) can be un- 
ambiguously determined. Because the three nucleons in 
this channel exist in a spin-symmetric state, and hence 



have an antisymmetric space-isospin wavefunction, the 
scattering in this state is completely determined by the 
long range part of the triplet s-wave NN interaction in 
the n-p channel; i.e. by n-p scattering and the properties 
of the deuteron. Furthermore, the use of multi-energy 
phase shift analysis, such as that conducted by Black et 
al. in the p-d system 20], could allow one to extrapolate 
to the very low energy quartet phase shifts on the ba- 
sis of data at higher energies. This is possible because, 
unlike the 2 Si/ 2 effective range function, the 4 S3/ 2 ef- 
fective range function is free of singularities. Therefore, 
one can independently predict the zero energy quartet 
s-wave scattering length 4 a n( j both from fundamentally 
sound theory and from the totality of the sub-breakup 
n-d database. 

A high accuracy determination of the 2 Si/ 2 scattering 
length in the n-d system is of interest for a number of 
reasons. The Pauli principle does not deter this channel 
from exploring the shorter-range components of the inter- 
action, where 3N forces should appear, and the Coulomb 
interaction is not present to complicate experiment or 
analysis. The n-d system is also in principle sensitive to 
charge symmetry breaking (CSB) effects. Due to charge 
symmetry breaking, the n-p force and the n-n force are 
different. CSB is explicitly included in modern NN po- 
tentials such as the Argonne AV18 potential. The accu- 
racy in this observable that can be achieved in modern 
calculations is sufficient to be sensitive to CSB. 

It was discovered in the early days of work on nu- 
clear three-body systems that there is an empirical rela- 
tion between the calculated doublet n-d scattering length 
and the binding energy of the tri-nucleon system. These 
so-called "Phillips lines" |2l| have been shown to be 
strictly linear for the n-p system, (in which case it 
is the triplet scattering length that is linked with the 
deuteron binding energy) and they are approximately 
linear for n-d scattering. When Phillips lines are plot- 
ted for the n-d calculations, the physical triton binding 
energy (8.48 MeV) intersects the currently accepted dou- 
blet value 2 a n d = 0.65 fm '22], in excellent agreement 
with the measurement of Dilg et al. of (0.65 ± 0.04) fm 
|23|. By contrast, a charge-symmetric R-matrix analy- 
sis of n-d data conducted by Hale 0] at energies up to 
the three-body breakup threshold gives the preliminary 
value 2 a n d = 0.41 fm, even though the Dilg et al. mea- 
surements were the only ones included in the n-d data 
set. These results, although preliminary, call into ques- 
tion the experimental values of the n-d scattering lengths, 
particularly in the doublet spin state. If the value of 2 a n d 
really is significantly less than 0.65 fm, then the presum- 
ably well-understood phenomenology of the Phillips line 
is thrown into disarray. 

The current values of the n-d scattering lengths were 
determined by measurements performed almost 30 years 
ago using a neutron gravity reflectometer |2^]. These 
measurements were not performed using pure D 2 gas 
samples because the experimental technique required 
samples in liquid or solid form: they were performed 
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on D2O, SiC>2 and Si from which the scattering length 
for n-d was deduced. This report describes the first 
high-precision measurement of the n-d coherent scatter- 
ing length performed using a pure sample. In addition, 
we chose to pursue a precision measurement of the n-p co- 
herent scattering length. The measurement of b np should 
yield a result that is consistent with the world average 
of other n-p measurements. If it doesn't then we may 
have an unaccounted for systematic effect in our mea- 
surement. Therefore the measurement of b np provides a 
secondary check on possible unknown systematic errors 
of our method. 

We use neutron interferometry to measure the coher- 
ent scattering lengths of these gases. Interferometry mea- 
sures a phase shift which is simply related to the coherent 
scattering length of the nuclei and is therefore sensitive 
to a particular linear combination of the scattering am- 
plitudes in the two channels. From an experimental point 
of view, it is possible to measure this linear combination 
to much higher precision than one can determine the in- 
dividual scattering lengths. For example, the coherent 
scattering length of silicon has been measured using neu- 
tron interferometry to the incredible accuracy of 0.005 %, 
surely one of the most accurate measurements of a nu- 
clear scattering amplitude ever performed|24|. The ex- 
periment requires neither polarized beams nor polarized 
targets and is sensitive to the amplitude as opposed to 
the square of the amplitude that enters into cross section 
measurements. 

The remainder of this paper is organized as follows. 
Section ITT1 outlines results from the theory of neutron op- 
tics that are relevant to the interpretation of the mea- 
surements. Sections IIII Al and IIII Bl describe the neutron 
interferometry technique and measurement. Section l!!! CI 
covers the gas cell design, which is also chosen to mini- 
mize the sample cell phase shift. Section fill Dl describes 
the procedure by which the target cell was aligned on 
the interferometer to make the measurement insensitive 
to the large phase shifts introduced by the gas sample 
cell. Section 1111 El presents the phase shift data, uncor- 
rected at this point for various slow changes in gas den- 
sity. Section UlI Fl describes the measurements which de- 
termine the gas density using the ideal gas law with virial 
coefficient corrections. Section UlI Gl shows the measure- 
ments which establish the gas composition of the deu- 
terium sample, which has non-negligible HD and H2 con- 
taminants. Section fill HI describes the sample thickness 
measurements with and without pressurized gas in the 
cell and presents the last time-dependent correction to 
the gas density, which removes all time dependence to 
the phase shift data. Section fill II outlines the measure- 
ment of the neutron wavelength and describes the mea- 
surement that establishes its stability over time. Sec- 
tion [W] relates our measurement of the coherent scatter- 
ing lengths of H2 and D2 molecules to the nuclear scat- 
tering lengths, outlines the correction to neutron optical 
theory due to the breakdown of the impulse approxima- 
tion, presents our final values for the scattering lengths, 



and compares them with past measurements. Section fVl 
describes the comparison of the world's data on the n-d 
coherent scattering length to existing theoretical calcu- 
lations and uses the recently-available high precision cal- 
culations of the quartet n-d scattering length in second 
generation NN potential models to determine the dou- 
blet n-d scattering length. It also places the measure- 
ments in the n-d system into the wider context of p-d 
measurements. Finally section I VII summarizes our re- 
sults and their significance and discusses the possibilities 
for further improvements in the precision of coherent n-p 
and n-d scattering length measurements as well as the 
possibilities for high-precision coherent scattering length 
measurements in n- 3 He, n- 4 He, and n- 3 H systems. Ap- 
pendixmdescribes in more detail the small corrections to 
the theory derived by Nowak, to which our experiment 
is sensitive. Appendix \2\ discusses the procedure used 
to determine the world average using published values of 
the scattering lengths. 



II. NEUTRON OPTICS THEORY 

In this section we give a brief review of the relevant 
results from the theory of neutron optics that are needed 
to understand precisely how the phase shift measured in 
the interferometer is related to the neutron scattering 
lengths of interest. For a more detailed treatment see 
Sears H. 

Neutron optics is based on the existence of the "co- 
herent wave" which is the coordinate representation of 
the coherent state formed by the incident wave and the 
forward scattered wave in a scattering medium. It is 
determined by the solution of a one-body Schrodinger 
equation 



2m 



A + v(r) 



ip(r) = Eip(r) 



(4) 



where ip(r) is the coherent wave and v(r) is the optical 
potential of the medium. The coherent wave satisfies the 
Lippmann-Schwinger equation 



if)(r) = |fc) + gv{r)ip(r) 



(5) 



where \k) is the incident wave, g is the one-body Green's 
function for nonrelativistic motion of a neutron and v(r) 
is the optical potential. The optical potential is related 
to the one body t matrix by 



t — v(r) + tgv{r) 



(0) 



and this combination forms the usual coupled system of 
equations of nonrelativistic scattering theory. Given a 
form for the t matrix one can determine the optical po- 
tential and then solve the one-body Schrodinger equation 
for the coherent wave. 

One must make an approximation for the t matrix of 
the neutron in a medium of scatterers. The usual approx- 
imation is essentially the Born approximation in which 
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v = t and 



(7) 



Finally one must approximate the one-body t matrix t;. 
Using the impulse approximation for scattering, one gets 



ti = (2vrft 2 /m)^6 i( 5(r-R i ). 



(8) 



Here, I denotes the elemental species, 6/ is the coherent 
scattering length for element I, r is a random spatial co- 
ordinate and R; defines the coordinate of each atom the 
neutron can scatter from. From Eq. IHJwe then arrive at 
an expression for the optical potential 



v opt (r) = (2TTh 2 /m)J2Nib l: 
I 



(9) 



where Ni is the number density of scatterers. 

The effect of the optical potential on the beam for a 
nonabsorbing uniform medium is to slow down the neu- 
trons as they encounter the potential step due to the 
matter, thereby decreasing the neutron wave vector, K, 
within the medium. The neutron index of refraction is 
defined by this relative change in the magnitude of the 
wavevector n — K/k. Conservation of energy at the 
boundary determines the relation to the optical poten- 
tial 



1 - v /E. 



(10) 



In general the scattering amplitude is complex to ac- 
count for incoherent scattering and absorption of the 
wave amplitude. The imaginary part of the scattering 
amplitude is related to the total reaction cross section by 
the optical theorem and leads to the more accurate form 
of the index of refraction 

n = n r + iri{ 



1 - (m\ 2 /2n) \jb 2 — (<7i/2X) 2 (11) 



+ iN m (A/4tt) 



Here, n r (rii) is the real (imaginary) part of the index 
of refraction, A is the neutron deBroglie wavelength, b is 
the bound coherent scattering length, and oi is the total 
reaction cross section (scattering plus absorption). For 
typical neutron-nucleus potentials, the real part of n is 
near unity, (n r — 1) w 10~ 5 . In less accurate treatments, 
the second term in the square root of the real part is 
omitted. However, its presence is required by the optical 
theorem and is included here for completeness. 

In neutron optics experiments with unpolarized neu- 
trons and unpolarized samples the coherent scattering 
length b is the sum of the scattering lengths in both scat- 
tering channels weighted by the number of spin states in 



each channel. From a quantum mechanical point of view, 
this is the total amplitude for a neutron to scatter with- 
out a change in the internal state of the target. In the 
forward direction only this amplitude can interfere with 
the unscattcrcd incident wave. In atomic deuterium with 
nuclear spin 5=1, for example, the coherent scattering 
length is the weighted sum of the doublet (S — 1/2) and 
quartet (S = 3/2) states 

]\/[ _i_ jji 

bud = [(1/3) 2 a nd + (2/3) 4 a» rf ] , (12) 

and for atomic hydrogen it is the weighted sum of the 
corresponding singlet and triplet states 



bn P = [(1/4) l a np + (3/4) 3 a np ] 



(13) 



Here the superscripts label the number of spin states for 
each scattering channel. In the next section the relation 
between the phase shift of the coherent wave in the in- 
terferometer and the bound coherent scattering length is 
derived. 



III. EXPERIMENTAL PROCEDURE AND 
RESULTS 

A. Apparatus and basic ideas of neutron 
interferometry 

Scattering length measurements were performed at the 
National Institute of Standards and Technology (NIST) 
Center for Neutron Research (NCNR) Interferometer and 
Optics Facility 261 . This facility, situated on cold neu- 
tron guide seven (NG7), consists of a perfect crystal 
silicon neutron interferometer with high phase contrast 
(80 %) and long-term phase stability (kj 5 °C per day). 
The single crystal interferometer is schematically illus- 
trated in Fig. ^ A monochromatic cold neutron beam 
(E = 11.1 MeV, A = 0.271 nm, AA/A < 0.5 %) 
is diffracted into the facility using a pair of pyrolytic 
graphite (PG) crystals operated in a non-dispersive dou- 
ble crystal mode using the (002) reflection. This beam is 
collimated and vertically focused onto the interferometer. 
The neutron beam is coherently divided in the splitter 
slab near point A by Bragg diffraction into two beams 
that travel along paths I and 77. These beams are again 
split in the mirror slab near points B and C by Bragg 
diffraction. Two of the wave fields are coherently recom- 
bined to interfere in the third Si mixer slab near point 
D. The perfect single crystal nature of the device ensures 
the required alignment precision which allows the Bragg 
diffraction condition to be met by all three silicon slabs 
simultaneously. In addition, the narrow phase space ac- 
ceptance of the device is an efficient filter for all incoher- 
ent neutron interactions in the silicon, phase shifter, and 
samples. 

The geometry traced out by this kind of interferome- 
ter is commonly referred to as Mach-Zehnder by analogy 
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FIG. 1: A schematic view of the Si perfect crystal neutron in- 
terferometer. Parameters associated with the neutron optics 
are discussed in the text. 



with similar optical interferometers. Here, however, the 
splitter, mirror and mixer crystals have been machined 
from a monolithic Si ingot in which the lower backbone 
maintains the perfect lattice registry of the original in- 
got. This avoids the mechanical difficulties associated 
with permanently aligning separated crystals to toler- 
ances much less than the angular acceptance of a per- 
fect crystal (5 /irad). The interferometer crystal geome- 
try is denoted as a Laue-Laue-Laue (LLL) crystal since 
there are three transmission Bragg reflecting crystals in 
the Laue geometry. For a detailed description of neutron 
interferometry techniques and experiments, refer to the 
book by Rauch and Werner |27| . 

The relative intensities recorded for the two exit 
beams, called the O-beam and the H-beam here, depend 
upon the phase difference of the neutron waves traversing 
path 77 relative to path I. The phase shift that we mea- 
sure is due to the phase difference of the neutron wave 
traversing path I with the sample in this path relative 
to its phase with the sample removed. For the samples 
used in these measurements, the attenuation cross sec- 
tion does not significantly influence the real part of the 
refractive index. To very high accuracy, the phase shift 
due to the sample can be written as, 



A(j) = (fir - l)kD eff = - WhDe 



(14) 



where D e g is the effective neutron optical thickness of the 
sample medium along the direction of wave propagation. 
Thus, a measurement of the coherent scattering length b 
to 0.05 % absolute accuracy demands an absolute mea- 
surement of the following quantities at the 0.02 % level: 



(1) the neutron optical phase shift of the gas, A<f>, 

(2) the atom density N, 

(3) the sample thickness D e g, and 

(4) the wavelength A. 

In addition, it is necessary to verify the purity of the gas 
at the same level of accuracy. Sections |lII Bl through IIII II 
discuss these measurements in detail. 



B. Measurement of the molecular H & D gas phase 

shift 



The phase shift, A</>, is measured by a secondary sam- 
pling method in which a flat phase shifter plate (denoted 
as the quartz phase flag in Fig. ^) is positioned between 
the splitter and mixer blades to intercept both neutron 
paths. This phase flag has an optical thickness described 
by 



AD flag = d[ 1 j 9 



m = 



D fiag = 2D f la9 £(A5), (15a) 
sin(0 B )sin(A<$) 



cos 2 (0 B ) -sin 2 (AS) 



Here Do is the thickness of this quartz flag (1.5 mm), 8b 
is the Bragg angle of the neutron beam reflecting from 
interferometer Si (111) lattice planes, and AS = (S— 5q) is 
the rotation angle offset of the phase flag in the horizontal 
plane. 

The intensities of the O and H beams (see Fig. ^| as 
a function of the phase flag angle, AS, are referred to 
as interferograms and can be described by the following 
relations: 



Io(S) 



RRTiPo + Te- iAm RRi>o 
A +B cos (C£($) + A#,) 



(16) 



and 



Ih(S) = 



TRTipo + Re-^^RRipo 



(17) 



.4 



H 



Bcos(C£(<5) + A<^+7r). 



Here R or T is the Si crystal reflection or transmission 
coefficient, and the parameters Aq, Ah, B, C, and A(f>' 
are extracted from fits to the data. The value of A<t>' 
and its corresponding uncertainty is used to determine 
the phase difference between the two interfering beams. 

This experiment measures the phase shift when a gas 
sample is placed into neutron path I (see Fig- EJ) • The 
gas and sample housing (cell) contribute separate phase 
shifts, A0q — > A(j) gas + A<f> cc u + A(j>Q. The effective thick- 
ness of the sample depends upon misalignment according 
to, 



D eff (e,j) = (A,)/[cos(Ae)cos(A 7 )], 



(18) 



where Dq is the true thickness and Ae = (e — eo) and 
A7 = (7 — 70) correspond to the horizontal rotation and 
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FIG. 2: Schematic view of the interferometer with gas cell and 
quartz alignment flag. The quartz alignment flag is shown in 
the out position and the centerline denotes an end on view of 
the kinematic mounting plane. The position of the tempera- 
ture probes (labelled Temp A and Temp B) are shown on the 
gas cell. 



vertical tilt of the sample. The condition e — £o and 
7 = 7o occurs when the beam is incident normally on 
the sample. In this experiment A(f> gas + A(f> ce u + A<fo is 
measured when the beam is incident normally. In this 
orientation the phase shift is insensitive in first order 
to the misalignment angles Ae and A7. However the 
phase shift due to the empty aluminum cell is two or- 
ders of magnitude larger than the gas. If the cell was 
present in only one of the interferometer beams, a small 
unknown misalignment would lead to a large systematic 
uncertainty in the background phase shift. To minimize 
this systematic effect a unique design was chosen for the 
gas sample cell, as discussed in the next section. 



C. Gas cell design 

In order to minimize the phase shift due to the cell, 
the cell walls were designed to extend across both beam 
paths to produce compensating phase shifts for paths I 
and II. Mechanically this was achieved by machining 
two gas cells from a single block of aluminum. The design 
of the cell is shown in Fig. [21 When the cell is perfectly 
aligned the beams strike both cell compartments perpen- 
dicular to their surfaces. This ensures that the effective 
thicknesses along both paths are nearly equal. The phase 
shift due to the cell obtained from Eas. ll4l andll8lis. 



A(j> ceU (e,j) 



-XNb 



{D H - Dj) 
cos(Ae) cos(A7) 



where D] and Djj are the cell wall thicknesses along 
paths I and 77. From Eq. ED we see that the phase shift 
along path / is opposite that of path II so that the total 
phase shift is minimized when Dj and Djj are equal or 
nearly equal, as is the case. In this experiment the total 
phase shift of the cell was first measured before filling 
with D 2 gas to be (A<j) ceU ) D2 = (2.4794 ± 0.0021) rad 
and (A(j> ceU ) H2 = (1.3788 ± 0.0021) rad. The two cell 
phase shift measurements differed by a slight amount due 
to thin (?« 10 /im thick) film of thermal grease that 
was present during the the D 2 gas measurement, but 
was removed (with acetone) prior to the H2 measure- 
ment. These cell phase shift values are considerably lower 
than 670 rad, which would have been seen with the cell 
in only one beam. The relative uncertainty introduced 
by a 17 mrad (1°) misalignment of the cell is 0.001 %. 
This is an order of magnitude below the uncertainty goal 
of this experiment. In practice, it is possible to align 
the cell in the beam such that both Ae < 2 mrad and 
A7 < 2 mrad so that the systematic alignment uncer- 
tainty is completely negligible. 



D. Alignment of the cell 

Although the experiment is insensitive in first order to 
both of the missets Ae and A7, it is still necessary to 
actually measure these values. The fact that the cell was 
designed to be extremely insensitive to changes in these 
angles introduces some experimental difficulties in actu- 
ally performing these measurements. Aligning the cell 
experimentally requires measurement of the cell phase 
shift as a function of both Ae and A7. To solve this 
problem a kinematic mount was designed to allow the 
cell to maintain the previous alignment Ae and A7 rela- 
tive to the beam after being removed and replaced. The 
alignment phase shift measurement was performed us- 
ing a quartz alignment phase plate that only crossed one 
beam in order to produce a noticeable phase shift. (See 
Fig- El) This optically flat, 1.5 mm thick quartz plate 
was mounted on the same flat surface as the cell. The 
procedure that was used to align the quartz is much the 
same as the alignment procedure described in Ref. [24| 
and is described below. However, here the sample is ro- 
tated by 7t/2 (90°) relative to the non-dispersive phase 
shift position discussed in p4| . 

This method of alignment required that the quartz 
alignment sample be parallel translated from path / to 
path II. Upon translation the horizontal misset angle 
changes sign (Ae — » — Ae) yielding a difference between 
the two phase shift measurements of 



e 7r/2 (e, 7 ) = A#e,7)-A#-e,7) 
XNbD 



(20) 



(19) 



cos(A7) cos(#b) 

x {2 + (Ae) 2 [l + 2tan 2 ((9 B )]}. 
Note that when A7 is held fixed, Eq. [201 is a quadratic 
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FIG. 3: Measurement of the difference phase Q^ii (Eq. I2()|l 
to obtain the minimum of the rotation alignment Ae and tilt 
alignment A7 as described in the text. 



function of Ae, which means that the center can be accu- 
rately determined in order to minimize the misset angle 
Ae. Similarly Ea. 1201 is a quadratic function of A7 when 
A7 << 1 and Ae is held fixed. Independently varying 
the tilt for a fixed value of Ae allows the experimental 
determination of the optimum tilt position for A7 = 0. 
The experimental measurements of these two alignment 
parabolas are plotted in Fig EI All phase shift data were 
taken with Ae and A7 at their minimum values. 



E. Data Collection 

Once the cell mount was aligned using the above tech- 
nique, the phase plate was replaced by the gas cell on 
its kinematic mount. The gas chamber on path / was 
then filled with sample gas, while the compensation cell 
oriented symmetrically on path II was evacuated. Inter- 
ferograms with the cell first in the "cell in" position and 
then translated to the "cell out" position (see Fig.^J) were 
collected in order to determine the phase shift due solely 
to the gas and the small difference in thickness between 
the two chamber walls of the cell (discussed earlier in 
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FIG. 4: A typical pair of interferograms denoting the change 
in intensity as the phase flag angle S is varied (see Fig. . 
Data are shown for both the cell filled with D2 gas in the 
beam and when the cell is removed from the beam path. 



section UlI C|) . Each set of interferograms (see Fig.0} re- 
quired approximately 42 min to obtain a measurement of 
A</> gas + A(j) ce u with a relative uncertainty of 0.3 %. The 
two interferograms in Fig.0]were fitted to Eas.HfilandH7l 
to obtain the time dependent phase shifts, A0o, shown in 
Fig- El The total phase shifts of the cell, the gas, and the 
time dependent offset phase shifts, A</>g OS +A</> ce ;;+A<^0) 
are plotted in Fig.EJ Since each measurement of A^o was 
performed within 42 min of the previous one, the time 
dependence of the phase shift could be directly measured 
and removed from the total phase shift. The time depen- 
dence of the empty interferometer phase shift, discussed 
later in section ITlI II is believed to be due to slight tem- 
perature fluctuations that cause small geometric shifts 
between the two paths of the interferometer. The phase 
shift values with the time dependence removed are plot- 
ted in Fig. [7| The total amount of data taken, 353 runs 
for D2 and 358 runs for H2, was based on a statistical 
uncertainty target of 0.02 %. 



F. Measuring the atom density N 

To obtain a value for the coherent scattering length 
from the phase shift data plotted in Fig. [7| using Eg. ITU 
several additional measurements must be performed and 
a few time-dependent density changes must be corrected 
for. First, it is necessary to determine the atom density, 
N, of the gas for each phase shift measurement. To do 
this we employ the following form of the virial equation 
for a real diatomic gas given by, 



N = 



2P 



k B T{\ + B P (T)P + C P (T)P 2 



(21) 



where P is the pressure, kg is Boltzmann's constant 
(value taken from [28|). T is the temperature, Bp is 
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FIG. 5: Empty interferometer (cell removed) phase shift, 
A</>o, plotted for each D2 (upper) and H2 (lower) run. 



the second pressure virial coefficient and Cp is the third 
pressure virial coefficient. (Note that the extra factor 
of two in Eq. [21 is included to account for the fact that 
each molecule contains two atoms.) The virial coefficients 
for hydrogen and deuterium have been measured |29l l30| 
with sufficient accuracy to determine N with a relative 
uncertainty of 0.001 % (for example: for deuterium Bp = 
5.616 x 10~ 4 bar" 1 , C P = 2.58 x 10~ 7 bar" 2 and for hy- 
drogen Bp = 5.796 xlO" 4 bar" 1 , C P = 2.36x 10~ 7 bar~ 2 
at T=297.15 K). The correction to N due to B P is about 
0.7 %, which is a significant contribution. The correc- 
tion due to Cp is about 0.004 % and is included in the 
calculation of N for completeness. 

To determine N using Eq. both the absolute tem- 
perature and the absolute pressure are needed. The ab- 
solute temperature of the gas was measured each time 
that the gas pressure was measured. To determine the 
absolute temperature, two calibrated 100 f2 platinum re- 
sistance thermometers were placed on the gas cell at the 
positions labelled Temp A and Temp B in Fig. [21 Al- 
though these thermometers have a precision of 0.002 %, 
the electronic readout accuracy was limited to 0.023 % 
at 300 K. These thermometers were calibrated at NIST 
by measuring the triple point temperature of pure H2O. 
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FIG. 6: Phase shift for the D2 (upper) and H2 (lower) filled 
gas cells as a function of measurement time. 

Before absolute calibration, probe A reported the triple 
point temperature as (273.057 ± 0.070) K and probe B 
reported (273.046 ±0.070) K, consistent within the abso- 
lute accuracy. Since the triple point of pure H2O is the 
operational definition of 273.15 K absolute temperature, 
a shift of (0.10 ± 0.07) K was applied to the temperature 
scales of both thermometers. The small differential non- 
linearity of the platinum thermometers allows us to apply 
the same shift at the 300 K temperature of the measure- 
ments. A plot of the temperature variation of the two 
probes using the original probe calibration curves cor- 
rected for the absolute calibration using the triple point 
is shown in Fig. El Note that for both gases the tempera- 
ture of the two cells was the same to within the absolute 
accuracy of the thermometry indicating that there is no 
thermal gradient across the cell. Since the thermal con- 
ductivity of the aluminum cell is high and there are no 
heat sources or sinks other than the surrounding envi- 
ronment the temperature of the gas can be inferred from 
the temperature of the aluminum cell by assuming ther- 
modynamic equilibrium. The data show there are slow 
temperature fluctuations: the amplitude of these fluctu- 
ations is consistent with the rated temperature stability 
of the interferometer enclosure. 

The pressure was measured using a Mensor DPG II 
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FIG. 7: Phase shift with the time-dependent initial phase shift 
subtracted off for (upper) D2 and (lower) H2 gas. The mean 
values and their first standard deviation are shown in the 
figure. One point in each figure is shown with the uncertainty 
estimated from Poisson counting statistics. 

model 15000 jU digital pressure gauge. This pressure 
gauge uses a silicon pressure transducer (SPT), which 
modifies the output from a precision micro-machined sil- 
icon wafer used as a pressure sensor. It was calibrated 
by the Mensor Corporation |31j using a dead weight test 
with NIST traceable standards. The dead weight test in- 
volves putting NIST traceable weights on a piston filled 
with dry nitrogen. The pressure exerted by the weight 
can be determined absolutely by knowing the masses and 
the cross sectional area of the cylindrical piston and com- 
pared with the output of the SPT. Calibration measure- 
ments with the dead weight established that the SPT 
readings were repeatable and linear functions of absolute 
pressure to better than 0.01 % |32(. 

This pressure gauge was directly coupled to the sam- 
ple cell through a high-pressure gas handling system 
schematically illustrated in Fig. |5J The top plate of 
the cell, which connects to the body through an indium 
o-ring seal, was connected to the gas handling system 
through two flexible 1 m long gas and vacuum lines made 
from HiP [3lJ thin stainless steel tubes 1.59 mm OD, and 
0.76 mm ID. The thin flexible lines also allowed the cell 



to be thermally insulated while allowing it to be trans- 
lated in and out of the interferometer without having to 
move the gas handling system. All connections in the 
gas handling system were made with all metal seals to 
minimize contamination of the sample. 

The pressure shown in Fig. 1101 decreases with time. 
This pressure change is not correlated with the temper- 
ature change of the cell nor is it due to an external leak 
into to the system from atmosphere. It is also not due 
to diffusion of the gas into aluminum: the known rate 
of diffusion of hydrogen into aluminum is much too slow 
|33|. This leak appeared only when one cell was evacu- 
ated and the other cell was pressurized. The conclusion 
drawn from this information is that the cell leaked into 
the evacuated cell from the filled cell through gas valve 
PV5 in Fig. |5J Although not an ideal condition, the pres- 
ence of such an internal leak does not compromise the 
purity of the gas. Also, the evacuated side was continu- 
ously pumped during the experiment so that there would 
not be enough gas in the evacuated section to cause a sys- 
tematic error in the phase shift measurement. The leak 
did change the atom density, N, at a rate of 0.002 % h -1 . 
However, this slow change was measured and can be cor- 
rected for easily in the data analysis. With data for both 
the temperature and pressure, the time-dependent atom 
density N is calculated from Eq. [21 with the results 
shown in Fig. ^2 



G. Gas Purity 

The D2 gas used in these measurements was "Research 
Grade" quality purchased from Air Products [3l]]. The 
quoted purity from the manufacturer was 99.99 %, but 
HD and H2 impurities were not included in this analy- 
sis. The H2 gas was "Ultra High Purity" purchased from 
Matheson Tri-Gas [3l| with quoted purity of 99.999 %. 

Measurement of the actual impurity species and con- 
centrations in the gas samples used is important to cor- 
rect the measured value of b at the 0.01 % level. To 
do this both mass spectroscopy and Raman spectroscopy 
were employed. The mass spectra allowed one to see the 
contaminating components in the samples as shown for 
example for D2 gas in Fig. ^| From the data, the main 
contaminants found in the D2 gas were HD at 0.3 % and 
D2O at a 0.02(1) % level. However, the relative ion- 
ization efficiencies of HD and D2 are not known for the 
mass spectrometer used in this measurement. Therefore, 
to measure the HD concentration to better than 5 % ac- 
curacy, Raman spectroscopy was used |34j . 

Raman spectroscopy measures the amount of light 
scattering from the rotational levels of HD and D2 , which 
are shifted due to the isotopic mass difference. Mea- 
surements were made using the apparatus described in 
Ref. [3^|. The spectrometer was first calibrated with a 
known 2 % HD sample (see Fig. 113b). Next, the D2 gas 
was introduced into the apparatus and rotational spectra 
were taken (Fig.ll3b.c). From this data the mole fraction 
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FIG. 8: The temperature of the cell measured with probes 
A and B (see text and Fig. |2J for each data run of the D2 
(upper) and H2 (lower) gas samples. The temperatures of the 
two probes agree within the absolute accuracy of the mea- 
surement. There is no evidence for a thermal gradient in the 
cell. 



of HD was determined to be x H d = (0.00301 ± 0.00013). 
The accuracy (4 %) of this measurement was limited by 
the sample pressure and in principle could be lowered to 
1 % with minor changes. 

With the Raman data for HD it is then possible to 
calibrate the mass spectrum data in order to determine 
the amount of other contaminants relative to D2. The 
mole fractions of D, H, and O were inferred to be xd — 
(0.99840 ± 0.00017), x H = (0.001500 ± 0.000065), and 
x = (0.000050 ±0.000016). The final expression for b D , 
corrected for impurities, was obtained using the following 
relation: 



bn = 



bgas — bHXH — boxo 
xd 



(22) 



where bn = 
0.004) fm H 



-3.7410 ± 0.0020) fm and b = (5.805 ± 
were used. 
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FIG. 9: Schematic of gas handling system used in experi- 
ments. During data collection valves PV3, PV4, PV6, and 
MV1 were open. 
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FIG. 10: Pressure of the D 2 and H 2 function of time. 

The pressure decreases slowly due to an internal leak into the 
evacuated side of the cell. A small correction for this effect is 
incorporated into the data analysis. 



H. Sample thickness 

Measurement of the sample thickness was performed 
using the NIST Precision Engineering Division Coordi- 
nate Measuring Machine (CMM) [37|]. This device is 
a temperature controlled coordinate measuring appara- 
tus capable of measuring macroscopic distances with an 
uncertainty of 250 nm within a 98 % confidence inter- 
val. This measurement showed that the thickness of the 
gas filled chamber of the cell was uniform at the 0.01 % 
level and that the cell was (1.0016 ± 0.0001) cm thick 
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FIG. 11: Atom density, N, of the D2 gas and H2 gas as a 
function of time. The atom density is calculated from the 
virial equation using the run-to-run pressure and temperature 
values Figs. 181 andllOl 



at an absolute temperature of (20.00 ± 0.05) °C. Varia- 
tion of the temperature resulted in a systematic change 
in the thickness at the 0.01 % level during the scatter- 
ing length measurement. Using the coefficient of linear 
expansion for aluminum a = 2.5 x 10~ 5 °C _1 38], this 
variation was taken into account in order to obtain the 
ratio Acj)/[N(1 + a AT)], where AT is the difference be- 
tween the temperature measured in the experiment and 
the temperature at which the cell dimensions were mea- 
sured by the NIST CMM. This ratio is plotted in Fig. HI 
from which a mean value and a standard error of the 
mean was calculated. 

Dimensional changes in the cell also occurred when the 
cell was pressurized to w 12.8 bar, which caused the walls 
to flex slightly outward. An upper limit for this effect 
was determined by both measurement using a Mitutoyo 
indicator [3lJ with 1 /an resolution and calculation using 
finite element analysis (see Fig. ll5ll . Both calculation and 
measurement confirmed that this flex resulted in a change 
in thickness at the center of the cell of less than 1 /im, 
which amounts to a systematic effect on the thickness of 
less than 0.01 %. 

The fact that the time dependence of the measured 
phase shifts shown in Fig. [21 has disappeared after 
all the known time-dependent effects on the gas den- 
sity are accounted for is a nontrivial result. The odds 
that other time-dependent effects not accounted for are 
present and conspire to cancel are very low, and there 
are no plausible physical mechanisms for other time- 
dependent effects on the phase shift beyond the ones 
mentioned so far (time-dependent wavelength shifts were 
eliminated as a possibility as described in the next sec- 
tion). Since the coherent scattering length being mea- 
sured is a time- independent constant, establishing the 
time-independence of the phase shift is a necessary con- 
dition for a valid measurement. 

In addition, the statistical spread of the measured 
phase shifts is consistent with that expected based on 
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FIG. 12: Low mass spectrum (upper) for D2 gas showing the 
relative peak heights of D + , HD + , and D^. The full mass 
spectrum (lower) shows the D3" peak height relative to all 
the higher mass contaminants. The ionization efficiencies of 
masses 18 thru 20 are about three times that of D2. Mass 19 
(fluorine) is a background due to outgassing of teflon in the 
mass spectrometer. 

Poisson statistics and the statistical accuracy of the indi- 
vidual phase shift measurements based on neutron count- 
ing statistics. This confirms that there are no unknown 
sources of noise present in the measurement and places 
further indirect constraints on possible interfering influ- 
ences on the data. 



I. Measuring A 

The wavelength was measured using an analyzer crys- 
tal in the H-beam of the interferometer shown in Fig. ^ 
Rotating the analyzer crystal through both the symmet- 
ric and the anti symmetric Bragg reflections allows the 
absolute Bragg angle, 9b, to be determined, and thereby 
the wavelength, using Bragg's law, A = 2 d sin^s), 
where d is the lattice plane spacing of the crystal. An- 
alyzer crystals of both pressed Si and pyrolytic graphite 
(PG) were used to determine 9b- Silicon is ideal since 
the lattice constant is known with an uncertainty much 
lower than our requirement of 0.01 The Si crystal 

was used at first along with the PG crystal to allow the 
PG lattice constant to be measured. The PG crystal was 
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FIG. 13: Raman spectroscopy of D2 gas. The upper plot is 
the calibration with a known HD ratio of 2 %. The middle and 
lower plots are the rotational spectrum of the D2 gas sample 
used in this experiment. This determined the calibration of 
HD to D2 peak in the mass spectrum data of contaminants. 
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mounted kincmatically to the interferometer setup to en- 
sure that the neutron beam was always sampled by the 
same part of the PG crystal. Therefore each measure- 
ment of the wavelength is correlated with the original 
measurement performed using silicon. 

The crystal was rotated about the vertical axis to lo- 
cate the two angles where the Bragg condition is satisfied. 
This condition is manifested as a dip in the transmit- 
ted beam measured in the 3 He detector. These rocking 
curves (intensity versus angle) are fitted to Gaussians to 
determine the centroids of the dips (Fig. I16f) . The wave- 
length, obtained from the 29 b angle was determined from 
the angular difference between the centroids. By per- 
forming measurements of the 29 b angle with both Si and 
PG with the same neutron wavelength the lattice con- 
stant of PG can be calibrated relative to that of Si. 

For the Si measurements it was necessary to include an 
additional linear term in the fit to account for a changing 
background due to a parasitic Bragg reflection of recipro- 
cal lattice points that lie near the Ewald sphere of reflec- 
tion. This effect does not appear in the data for PG since 
the lattice of PG is randomly oriented about the hexag- 
onal c-axis. (This random orientation effectively makes 
the reciprocal lattice of the PG crystal one-dimensional.) 

The measurement of the 29b angle is also sensitive to 
the relative tilt between the interferometer lattice planes 
and the analyzer lattice planes. In order to align the 
analyzer crystal, measurements of 29b were performed 



FIG. 14: Point-to-point gas phase shift divided by the point- 
to-point atom density, N, and the point-to-point correction 
for the length of the cell due to thermal expansion and the 
temperature difference between the neutron measurement and 
the length measurement. The mean values and their first 
standard deviation are shown in the figure. One point in 
each figure is shown with the uncertainty estimated from 
Poisson counting statistics. The final standard uncertainty 
(er/yN — 1) of the mean must be added in quadrature to the 
systematic uncertainties due to the pressure and temperature 
measurement. 



for various tilt angles of the analyzer lattice planes. The 
functional form of 29b plotted versus the tilt angle is 
a parabola (see Fig. I16f) with a minimum corresponding 
to the condition that the lattice planes of the analyzer 
crystal are perpendicular to the scattering plane of the 
interferometer. This method allows 9b to be determined 
with an uncertainty of 0.001 %. The resulting value for 
A, determined from the data shown in Fig. 1171 was found 
to be (0.271266±0.000012) nm for the D 2 gas experiment 
and (0.2713050 ± 0.0000085) nm for the H 2 gas experi- 
ment. These two numbers differ slightly at the 0.01 % 
level due to minor changes in the mechanical configura- 
tion of the apparatus between the experiments. These 
changes were due to the need to remove the cell between 
the measurements. 

Not only is it necessary to know the mean value of the 
wavelength of the beam passing through the interferome- 
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FIG. 15: The three-dimensional perspective of the gas cell modelled using finite element analysis is shown. Although the entire 
surface was mapped in the actual calculation, only the central and horizontal axes (where the maximum deflection occurs) are 
shown for a pressure of 13 bar. 



ter, but it is also necessary to ensure that this value does 
not change during the course of the experiment. Unlike 
the temperature and pressure measurement the wave- 
length was measured only twice during the experiment 
for practical reasons. To prove that the time-dependent 
phase drift described previously was not due to time- 
dependent wavelength changes a thick aluminum sample 
was placed in the beam which amplified the phase shift 
in the interferometer by a factor of 100, thereby ampli- 
fying any fluctuations in the phase due to possible wave- 
length changes by the same factor. The result of this test 
showed that as long as no mechanical changes in the ex- 
perimental configuration have been made, the long-term 
stability of the wavelength is better than 0.001 %. This 
means that, within limits of the statistical uncertainty, 
the phase drift is due to temperature fluctuations and 
not due to a drift of the wavelength. 



IV. RESULTS FOR b D AND b H 

All of the measurements required to calculate the co- 
herent scattering lengths in H2 and D2 gas are now de- 
termined. It remains to relate the coherent scattering 



lengths of the molecules, and bjj, to the coherent 
scattering lengths of the nuclei, b np , and b n d- 

The general expression for the low energy neutron scat- 
tering length of an atom away from nuclear resonances 
for unpolarized atoms and neutrons is 

b = b nuc + Z(b ne + b s )[l - f(q)} +b s + b P oi (23) 

where b nuc is the scattering amplitude due to the 
neutron-nucleus strong force, b ne is the neutron-electron 
scattering amplitude due to the internal charge distribu- 
tion of the neutron, b s is the Mott-Schwinger scattering 
due to the interaction of the magnetic moment of the 
neutron with the v x E magnetic field seen in the neu- 
tron rest frame from electric fields, b po i is the scattering 
amplitude due to the electric polarizability of the neutron 
in the intense electric field of the nucleus, and f(q) is the 
charge form factor (the Fourier transform of the electric 
charge distribution of the atom). We first note that the 
electromagnetic contribution to the scattering lengths of 
both hydrogen and deuterium from b ne and b s are exactly 
zero for forward scattering due to the neutrality of the 
atoms, which forces the charge form factor f(q) 1 as 
q — > 25]. In addition, the contribution from the elec- 
tric polarizability of the neutron for both H and D is less 
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FIG. 16: A typical set of rocking curves (intensity vs. angle) 
for the analyzer crystal at a fixed tilt angle (above). The 29 b 
angle is calculated from the difference between the center of 
the parallel Bragg curve and the anti-parallel curve. Below, 
the 29b angles are plotted as a function of the tilt angle of the 
analyzer . The minimum of the parabola corresponds to the 
correct analyzer tilt position, and therefore the correct Bragg 
angle 6b- 



than —0.000017 fm[40| and is therefore negligible. Thus 
this measurement is sensitive only to that part of the 
scattering length due to the strong interaction, which is 
precisely what can be calculated using theoretical models 
of the NN interaction. Although there are in principle ad- 
ditional contributions due to local field effects that arise 
from multiple scattering in the medium and its modi- 
fication of the amplitude of the incident waves on the 
scattering centers, these corrections are completely neg- 
ligible for gase s and only reach the 10~ 4 level in solids 
and liquids 25] . Therefore both our measurement, and 
the most accurate of the past measurements, measure 
the same physical quantity and therefore can be directly 
compared to each other as well as to NN theory. 

A summary of the parameters discussed in the pre- 
vious section are given in Table [I] Using these val- 
ues as input to Eq. H3 we obtain the average molecu- 
lar bound coherent scattering length per atom for deu- 
terium of bp = (6.6649 ± 0.0040) fm and for hydrogen of 
b H = (-3.7458 ± 0.0020) fm. 
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FIG. 17: Set of Bragg angle values for a set of measurements 
of the type shown in Fig. 1161 Note that the data is shown for 
the calibrated PG (002) analyzer crystal at the optimum tilt 
position. The value of the PG lattice parameter, c, calibrated 
against the Si-lattice parameter (a = 0.543101993 nm) is c = 
0.670982 nm. 



To compare this result with other measurements we 
must calculate an average value from the previous world 
values of the b np and b nc i bound nuclear coherent scat- 
tering lengths. The details of this calculation are left 
to Appendix an d the result of this evaluation shown 
in Fig. HU is b np = (-3.7410 ± 0.0010) fm and b nd = 
(6.6727 ± 0.0045) fm. Our result differs from the average 
of previous measurements by 2.3cr for H and 1.9a for D. 

This comparison is not yet complete since there re- 
mains one class of effects, which although not generally 
taken into account in past measurements, can in princi- 
ple be large enough to introduce a nontrivial correction 
in the usual relation between the scattering length of the 
molecule and the scattering lengths of the constituent 
atoms. In 1982, Nowak |41j revisited the approxima- 
tions inherent in the use of the Fermi pseudopotential 
to describe the refractive index, n, for neutrons, which 
had first been investigated by Lippman [42l | and was the 
first application of the Lippmann-Schwinger formulation 
of scattering theory 43] . Nowak came to the conclusion 
that the static local field corrections ^ were of the or- 
der of 10~ 4 , in agreement with Lax's earlier estimates 
[idl ]. However, he found that dynamic effects due to vir- 
tual excitations to the low lying states, within second 
order perturbation theory, make a relative contribution 
of order 10~ 3 to (1 — n). The effect is most pronounced 
for light target molecules such as H2 and presumably for 
this reason have rarely been taken into account in n-p 
coherent scattering length measurements, most of which 
involve measurements in which the hydrogen is bound 
in hydrocarbons. (These corrections are also mentioned 
in the n-p parahydrogen cross section measurements of 
Callerame |45j, which is one of the measurements used 
to determine the n-p scattering length.) The lowest 
lying excitations of the molecule are rotational states, 
E ro t = (S- 2 /2^)0)(j + 1), where / is the molecule's mo- 
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TABLE I: Parameters and relative uncertainties required to determine the scattering length. The coherent scattering length is 
calculated using the expression in Eq. 1141 the parameters in this table, and the correction for D2 gas composition in Eq. 1221 
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merit of inertia. For H2, (ti 2 /2I) — 7.56 meV; the vibra- 
tional energy levels are separated by 546 meV. Nowak 
carried out a detailed numerical calculation of this dy- 
namic correction, finding it to be a relative correction 
of 1.1 x 10~ 3 . A later estimate by Summerfield 
agrees with this conclusion. Nowak also calculated the 
correction to (1 — n) due to multiple scattering of the 
neutron within the H2 molecule. This leads to a relative 
correction of about 0.8 x 10 -3 . Both of these correc- 
tions are spin-dependent so that the para- and ortho- 
states of the H 2 molecule must be considered separately. 
The absolute correction to the scattering length calcu- 
lated from the sum of these two relative corrections is 
Ab Nowak = 1.9 x 1CT 3 x (-3.74 fm) = -0.007 fm, which 
is consistent with the difference between our neutron in- 
terferometry measurement of b np and the previous world 
average calculated using bn — b n d + Ab^ owa k- This pre- 
vious world average includes results from reflectometry 
measurements of hydrocarbon liquids. It is presumed 
that the measurements of b np for hydrogen bound in a 
hydrocarbon sample will not need to be corrected for 
multiple scattering or virtual excitations due to the in- 
verse molecule mass dependence of the effect. 

This measurement of b np by neutron interferometry on 
H2 gas may therefore be regarded as the first observation 
of the corrections to the index of refraction due to vir- 
tual excitations and multiple internal scattering within a 
molecule as predicted by Nowak 20 years ago. This con- 
clusion must be regarded as tentative at this point, since 
the calculation of Nowak was done in the long wavelength 
limit kR -^0(R = bond length of H 2 = 0.74611 nm, for 
D 2 Rq = 0.74164 nm j^)- We are currently performing 
a calculation for the conditions of our experiment, which 
correspond to /ci?o=1.73. 

If we apply the existing Nowak correction to our H2 
data, our result for the coherent n-p scattering length 



is b np = (-3.7384 ± 0.0020) fm. Taking this value 
we calculate the new world average value to be b np = 
(-3.7405 ± 0.0009) fm, (our result lies within la of the 
average). The Nowak correction is expected to be much 
larger for H2 than for D 2 due to the fact that there is a 
sign difference between the singlet and triplet scattering 
lengths of hydrogen, which results in an amplification 
of the correction to the scattering amplitude given by 
Eq. [213 The consequence of this sign difference is most 
prominently seen in the the incoherent scattering cross 
section of hydrogen, which as a consequence is dispro- 
portionately larger than the incoherent scattering cross 
section of deuterium. Therefore, we do not make a sim- 
ilar Nowak correction to the D2 data. With this uncor- 
rected value of b nc i we arrive at a new weighted average 
of b nd = (6.6683± 0.0030) fm (our result is within 0.9ct of 
the average value). Our results are consistent with past 
measurements conducted decades earlier by completely 
different methods. 



V. COMPARISON TO THEORY AND 
IMPLICATIONS FOR THREE-BODY FORCES 

We propose to show in this section that the physics 
impact of these coherent scattering length measurements 
in the n-d system are greater than has been recognized 
in the theoretical few body community. The reason is 
simple: rather than compare theoretical calculations to 
the coherent scattering lengths, which are measured to 
high precision, theorists in the field have instead com- 
pared their results to the scattering lengths in each of 
the two channels individually despite the fact that these 
are known with a precision that is lower by an order of 
magnitude. 

To illustrate the impact of comparing theoretical cal- 
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FIG. 18: Bound coherent neutron scattering lengths for n- 
d and n-p along with reported uncertainties. Our result for 
D2 is consistent with the current world average. Our result 
for H2, uncorrected for the multiple scattering and second- 
order perturbation theory corrections discussed in the text, 
is 2.3a away from the current world average. Including these 
theoretical corrections brings our value into closer agreement 
as shown. 



culation to the precision data on the coherent scatter- 
ing length, which has actually existed now for decades, 
we show several modern calculations of the scattering 
lengths in the p-d and n-d systems, some of which include 
a three-nucleon force. The dependence of the theoreti- 
cally calculated 2 Si/ 2 scattering length on the inclusion 
of a 3N force is clearly seen in Table [HI which gives the 
n-d and p-d doublet and quartet scattering lengths (in 
fm) calculated with a number of different potential mod- 
els. The results in boldface come from models in which 
the 3N force has been adjusted to replicate the trinucleon 
binding energy. 

In addition, the values of the coherent scattering 
lengths from Table [D] are plotted in Fig. ^| along with 
the world average value for b nc [ with la and 2a confi- 
dence bands. We note that, as expected, none of the 
theories which do not incorporate a 3N force of some 
sort come close to matching the n-d coherent scattering 
length. Of the potential models that include a 3N force 



TABLE II: Theoretical calculations of the p-d and n-d s- 
wave scattering lengths. In the last column we have cal- 
culated the coherent scattering lengths for the n-d system 
by forming the appropriate linear combination. We observe 
that none of the theories, with the possible exception of the 
AV14 potential with a 3N force, are in agreement with the 
precisely-known world average coherent n-d scattering length 
of (6.669 ± 0.003) fm. 
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of some tvpe, only the AV14 potential with the Brazil 
3N force[22j and the AV14 potential with the Tucson- 
Melbourne 3N force[5l] are in agreement with the data. 
The AV18 potential with the Urbana 3N force, the AV14 
potential with the Brazil 3N force, the Malfliet-Tjon- 
III potential, the Reid soft-core+Tuscon-Melbourne 3N 
force, the Reid soft-core adjusted to fit the triton bind- 
ing energy, and none of the ot her p otential models which 
include the Brazil 3N force ps. I5C| are in agreement with 
the data on the n-d coherent scattering length or, in cases 
where there is no separate calculation of the quartet n- 
d scattering length, in agreement with the doublet n-d 
scattering length derived below. This indicates that the 
precision with which the coherent n-d scattering length 
is known sets a tight constraint on NN potential models 
and, in our view, ought to be considered as automatically 
as is the value of the triton binding energy when such 
potentials are compared to the low energy observables. 
With the new result of the analysis of n-d scattering in 
chiral effective field theoryj^ including chiral 3N forces 
which shows that only two low energy parameters are 
required at NNLO order (surprisingly none are required 
at NLO|H3), the question of which observables to use to 
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TABLE III: Theoretical calculations of the quartet n-d s-wave 
scattering length. Table taken from Friar et al. 55] 
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fix these low energy constants becomes timely. We be- 
lieve that the best two low energy observables in the n-d 
system to determine these low energy constants are not 
the triton binding energy and the doublet n-d scatter- 
ing length but rather the triton binding energy and the 
coherent n-d scattering length. 

We can also use the coherent scattering length data in 
combination with theoretical calculations of the quartet 
scattering length to infer the doublet scattering length 
with significantly higher precision than the currently 
quoted value of 2 a n d = (0.65 ± 0.04(expt)) fm. As is 
well-known, the quartet scattering length is mainly sensi- 
tive to the well-known long-range components of the NN 
interaction due to pion exchange, which are fixed by mea- 
surements of the s-wave component of the deuteron wave 
function. So the results of NN model calculations should 
give the same answer for this channel to high accuracy 
independent of the details of the short-range components 
of the NN interaction where 3-body forces start to mani- 
fest themselves. Then this procedure can be used to more 
tightly constrain the short-range interactions, including 
3-body forces, which must be introduced to agree with 
the triton binding energy and to calculate the doublet 
scattering length. 

Table Hill shows calculations of the quartet n-d scatter- 
ing length using a new class of potentials which provide 
good fits to the NN database 55] . We observe that all of 
these results fit within a range A a n d = (6.346 ±0.007) fm. 
If we accept this average and range as a fair represen- 
tation of the precision with which modern NN poten- 
tials can calculate 4 aw, then we can combine this result 
with our measurement of the coherent scattering length 
to obtain the following value for the doublet scattering 
length: 2 a nd = (0.645±0.003(expt) ± 0.007(theory)) fm. 
Compared with the direct measurement 2 a n d = (0.65 ± 
0.04(expt)) fm, we see that this approach can improve on 
the precision of our knowledge of the doublet n-d scat- 
tering length by a factor of 4. 

At this point it would also be interesting to compare 
these new results with the scattering lengths in the p- 
d system. Unfortunately, there are at present no cur- 
rently accepted values for the low energy p-d scatter- 
ing lengths. The analysis of Black et al., did reproduce 
some of the theoretically predicted energy dependence of 
the p-d 2 Si/2 effective range function-a singularity near 
threshold-but the singularity was located at higher than 
expected energies and did not match the theoretical val- 
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FIG. 19: Theoretical calculations of the coherent scattering 
length compared with the experimental value measured here. 
The central dark band is the la confidence band and the 
lighter band is the 2a confidence band. Only one of the the- 
ories fall within the la band and only 2 fall within 2a. 

ues within quoted uncertainties |2(i| . None of the other 
experimentally determined values for this parameter were 
even close to the calculated values, and the effective range 
functions did not display the expected pole. For the p-d 
case, the conflict between theory and experiment persists 
whether or not 3N forces sufficient to correctly bind 3 He 
are added to the potential. 



VI. CONCLUSIONS AND FURTHER WORK 

We have performed high precision measurements of 
the coherent neutron scattering lengths of gas phase 
molecular hydrogen and deuterium using neutron in- 
terferometry. We find b H = (-3.7458 ± 0.0020) fm 
and b D = (6.6649 ± 0.0040) fm. Our result for H dif- 
fers from the world average of previous measurements, 
b np = (-3.7410 ± 0.0010) fm, by -0.0048 fm, which 
is accounted for by the Nowak correction. Our re- 
sult for the n-d coherent scattering length is in agree- 
ment with the world average of previous measurements, 
bnd = (6.6727 ± 0.0045) fm. We feel that the precision 
of these results has yet to be properly appreciated. We 
note that calculations of the doublet and quartet scat- 
tering lengths of the best potential models show that 
almost all known calculations are in disagreement with 
the precisely-measured linear combination correspond- 
ing to the coherent scattering length. Combining the 
world data on bo with the modern high-precision theo- 
retical calculations of the quartet n-d scattering lengths 
recently summarized by Friar et al., we deduce a more 
precise value for the doublet scattering length of 2 aw = 
(0.645 ± 0.003(expt) ± 0.007(theory)) fm. This value is 
a factor of 4 more precise than the previously accepted 
value of 2 a n d = (0.65±0.04(expt)) fm and is in agreement 
with the Argonne AV18 potential with a 3N force. 

We hope that this work will contribute to the exten- 
sive theoretical and experimental efforts now underway 
to understand the nuclear 3-body force. When possible, 
we urge theorists who calculate both scattering lengths 
in the n-d system to compare to the precisely known co- 
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herent scattering length in addition to the lower-precision 
values of the separate scattering lengths. In trying to use 
the n-d system to constrain possible forms of the 3-body 
force this procedure should be more sensitive. In partic- 
ular, we hope that the n-d coherent scattering length and 
the triton binding energy will be used to constrain the 
two NNLO low energy constants that model 3N forces in 
chiral effective field theory. 

It would of course be very useful to have another pre- 
cision measurement in the n-d system which is sensitive 
to a different linear combination of scattering lengths so 
that both scattering lengths could be extracted from ex- 
periment alone with high precision. Performing a inter- 
fcrometry measurement with a polarized D 2 target of suf- 
ficient density to operate stably in a neutron interferome- 
ter seems impractical at the moment due to the cryogenic 
and high magnetic field requirements needed to make 
polarized D2. Although polarized D2 gas targets using 
spin-exchange optical pumping are under development, 
the densities reached so far are still too low for a precise 
neutron interferometer experiment |59j| . One possibility is 
to perform a measurement of the pseudomagnetic preces- 
sion of polarized neutrons in a polarized D 2 target. Such 
a measurement determines a phase shift proportional to 
the difference in the scattering lengths. Recently such a 
high-precision measurement was performed with polar- 
ized 3 He at the ILL 60]. It would be very interesting if 
it were possible to adapt the elegant neutron spin pre- 
cession measurement technique used by Zimmer et al. in 
this measurement to the case of polarized D or H. 

Another natural question to ask is how the current 
measurements can be improved if more neutron flux with 
a neutron interferometer were to become available. The 
absolute measurement of the gas density can be improved 
by at least an order of magnitude before one needs to 
worry about the accuracy of the virial coefficients. A 
skew symmetric interferometer with 10 cm beam paths 
could be constructed that would allow the sample cell 
path to be increased by nearly an order of magnitude 
reducing the relative dimensional uncertainty by that 
amount as well. Designing the interferometer to work 
at a longer wavelength of 0.38 nm would also increase 
the sensitivity by a factor of 1.4. Measurement of the 
HD contamination of the gas in the case of D2 could be 
improved from 4 % relative uncertainty to 1 %. With 
all of these improvements taken into account the uncer- 
tainty can in principle be reduced by nearly an order of 
magnitude. 

Finally it is worth noting that the last decade has wit- 
nessed exciting advances in the accuracy of calculations 
in few-body nuclei. Therefore it is possible now to envi- 
sion the accurate calculation of low energy neutron scat- 
tering lengths for systems with A > 3. The next most 
interesting systems to measure are clearly 3 H, 3 He, and 
4 He, and as room temperature gases all the techniques 
described in this paper are applicable to them. In a fu- 
ture paper we will present results of a precision measure- 
ment of the n- 3 He coherent scattering length, which in 



combination with the recent measurement of the scatter- 
ing length difference in the same nucleus will be used to 
determine separately the scattering lengths in both chan- 
nels. This should improve the accuracy of our knowledge 
of low energy scattering observables in the 4-body system 
by an order of magnitude. It would also be possible to 
improve the current measurements of the coherent neu- 
tron scattering length of tritium. Again the motivation 
would need to come from theoretical developments in the 
4-body problem. 

Any NN potentials that correctly describe the two and 
three-body systems, even if a three-body force needs to 
be added, must predict correctly the n- 3 He scattering 
length with no adjustable parameters. A deviation would 
indicate either the presence of four-body forces (but these 
are believed to be even smaller than three-body forces on 
theoretical grounds) or a distortion of the NN interaction 
in the nuclear medium. Recent theoretical calculations 
exist for the n- 3 He system [6l|. 
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APPENDIX 

1. Corrections to the Impulse Approximation used 
in Neutron Optics and Consequences for High 
Precision Scattering Length Determinations 

We have chosen to describe in some detail the assump- 
tions of the usual theory of neutron optics in section [n] 
because we need to consider an effect which violates some 
of the approximations which undcrly this formalism. The 
approximation fails in two places: in Eqs. 0and|S| due 
to effects from atomic binding and multiple scattering 
respectively. 

To understand the effects of atomic binding, it is worth 
recalling the approximations which underly equation [S] 
Since the neutron-nucleus interaction is much stronger 
and much shorter range than the binding forces of the 
atoms in matter, it is reasonable to neglect the effects of 
chemical binding during the neutron-nucleus collision. In 
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addition, the short-range of the interaction means that 
the timescale of the collision is much shorter than the 
timescales associated with the motion of the atom in the 
potential well. For both of these reasons, the t operator 
is usually approximated by the t operator for a free atom. 
This is known as the impulse approximation in scattering 
theory. 

In this approximation, the neutron optical properties 
of a medium depend only on the coherent scattering 
length of the atoms and not at all on the details of 
the binding of the atoms. This independence from the 
medium has been exploited for decades in precision mea- 
surements of neutron coherent scattering lengths, which 
have typically been determined by measurements in mix- 
tures and compounds chosen for their experimental con- 
venience. As a matter of principle, however, at some 
point this approximation will no longer be valid and bind- 
ing effects will need to be taken into account. When this 
regime is reached, it will be necessary to know enough 
details about the atomic binding to calculate their ef- 
fects on the optical potential. It is our assertion that 
this regime has now been reached in our experiment. 

The calculation of corrections to the impulse approx- 
imation for neutron scattering lengths was plagued for 
decades by the singular nature of the delta function po- 
tential for the t matrix |42l Il2l 1621 l63l l64j . However in 
a series of papers |4l| l65l l66| manifestly divergence- free 
expressions for the first-order corrections to the impulse 
approximation were finally obtained which are consistent 
with the optical theorem and reduce in appropriate limits 
to previous results in the limit of static scatterers [63 . As 
an example, the case of a nucleus of mass M in an atom 
bound in a harmonic oscillator potential, the first-order 
correction to the impulse approximation was evaluated 
analytically 65] with the result 
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where / is the scattering amplitude, A = M/m is the 
ratio of the nuclear mass to the neutron mass, b is the 
bound scattering length, and a is the average amplitude 
of the oscillator in the ground state 
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with ui the oscillator frequency. The scale of the correc- 
tion is set by the ratio of the scattering length (a few 
fm) to the vibration amplitude of the atom (0.1 nm or 
so), which is small and typically neglected in neutron 
optics theory. However, since the size of this correction 
is a decreasing function of A, it is largest for H and D. 
In the particular case of a bound proton the scattering 
amplitude becomes 
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where 63 and b\ are the triplet and singlet n-p scat- 
tering lengths and the numerical factor of 0.578 comes 



from evaluating the integral. The factor multiplying the 
expression is the usual coherent scattering length for 
hydrogen. For the particular case of hydrogen, where 
63 = 5.42 fm and 61 = —23.75 fm, the size of the correc- 
tion term is anomalously large from the small denomina- 
tor, 1.0 x 10~ 3 for a typical atomic oscillation amplitude 
of 0.05 nm. 

However there is another effect of the same order of 
magnitude that must be taken into account associated 
with intraparticle multiple scattering. In Eq.[7|the t ma- 
trix of the bound system of N scatterers is expressed as 
the sum of the (impulse approximation) one-body t ma- 
trix. But this is known in exact treatments of scattering 
theory to be an approximation!^. The next order of 
approximation for the t matrix of the system is 



* = !>+ E tlGt r 



(27) 
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where G is the Green's function. Therefore the t operator 
for the system to second order is not the sum of the t 
operators for the individual particles due to interparticle 
multiple scattering. So to calculate the full t operator of 
the system, which is what is required to obtain the optical 
potential in Eq. 0] one must take multiple scattering also 
into account. Then finally the optical potential, which 
is now a function of the neutron momentum, must be 
solved from Eq. [B] and inserted into Eq. 1101 

Nowak^l] performed this calculation to obtain the 
modified expression for the index of refraction 



nf = 1 - 



E 



(28) 



and two correction terms at second order: one from the 
binding potential 



+00 



, ( , 2) = A dt d 3 ae^ k2 -^ 



(29) 



1 



,2 e i(q-k)Ri(t) e i(k-q)R,(0) 



and the other from multiple scattering 
v$ = A I dt I d\e^ k2 -^ (30) 



1,1', tyi' 

6;el (q-k)R i (t) \ U ^(k-qJR,, (0) 



and 



.4 



-tp 



{2-k) 6 h 



2ttH 2 



(31) 



In these equations p is the number density of scatterers, k 
and q are the wave vectors of the incident and in- medium 
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neutrons (very close to identical for cold neutrons), bi is 
the (spin dependent) scattering length operator of atom 
I, and the averaging () is the usual trace over spins and in- 
ternal states of the scattering system. These expressions 
are related to (but not identical with) the well-known dy- 
namic structure factors defined in the theory of thermal 
neutron scattering|69] and make it clear that the second 
order approximation to the optical potential is a func- 
tion of the dynamics and correlations of the scattering 
medium. 

From this expression we can understand why the mul- 
tiple scattering correction to the optical potential is es- 
pecially large for the case of H and D gas molecules and 
much smaller for hydrogen atoms embedded in the large 
polyatomic molecules that have been used for the most 
precise n-p and n-d coherent scattering length measure- 
ments in past work in liquid. First of all, the multiple 
scattering term vanishes if the nuclei in the poly- 
atomic molecule possess uncorrelated nuclear spin direc- 
tions and uncorrelated relative motions, which is true 
to an excellent approximation. From a physical point of 
view, multiple scattering from uncorrelated nuclei cannot 
contribute to coherent scattering. By contrast, as is well- 
known, the spins of hydrogen and deuterium molecules 
are tightly correlated as a result of the identical parti- 
cle constraints, and the relative motions are perfectly 
correlated for isolated molecules by the conservation of 
momentum. Therefore the multiple scattering correction 
will be much larger for H and D molecules than for larger 
polyatomic molecules. 

In addition, one can see that the size of this effect will 
be larger for hydrogen than for deuterium. In addition 
to the effect of the larger mass, deuterium has a smaller 
ratio of incoherent to coherent scattering, and as one can 
see from Eq. the effect is proportional to this ratio, 
which is anomalously large for hydrogen. 

The expression for the dynamic structure factor in the 
corrections calculated by Nowak can be evaluated exactly 
for hydrogen and deuterium gas in thermal equilibrium in 
the rigid rotor approximation. The calculation draws 
on the work of Young and KoppeljjJl on the dynamic 
structure factor for neutron scattering from hydrogen and 
deuterium molecules. The rigid rotor approximation is 
known to be an excellent approximation for slow neu- 
tron scattering in hydrogen and deuterium. Since the 
neutron energy is too low to excite the first vibrational 
level only transitions to rotational states need to be taken 
into account. These detailed calculations are presently in 
progress. 



2. Calculating the world scattering length average 



ments of the coherent scattering length of h ydr ogen 
MMMMMMMM and deuterium]!! E3, 

We excluded all mea- 
surements that were not published in a refereed jour- 
nal and all measurements that were later retracted for 
hydrogen[Hli3 and for deuteriumjlliliillJillli 
with the exception of Ref. which updates a previ- 
ous value reported in Ref. |87j for isotopic purity. Al- 
though the result of Bartolini et aZ.[9(J was never for- 
mally withdrawn it was excluded here due to the fact 
that the result was > 10c from the world average. The 
value quoted by Bartolini et al. 90] was due to a re- 
analysis of the H2O content of the sample used 4 years 
prior in Ref. (s^. This large discrepancy was also dis- 
cussed by Dilg et al.%2^, and Nistler who found this 
result to be inconsistent. In addition, there were two 
measurements of the ratio of the bound coherent scat- 
tering lengths of hydrogen and carbon whose originally 
reported uncertainties could be lowered by a factor of 2 
as a result of subsequent precision measurements of the 
bound coherent scattering length of carbon, which is now 
known to be b c = (6.6484 ± 0.0013) fm'ST]. Including 
this new value for be with the measurement of Dickinson 
et d.[z3 of b c /b H = (-1.775 ± 0.004) fm gives b H = 
(—3.748 ± 0.009) fm for which permission was obtained 
from one of the coauthors |94| to modify the original pub- 
lished result (original value: bn = (—3.74 ± 0.02) fm). 
However permission was not obtained to modify the mea- 
surement of Burgy et q/.|72|. Without permission to 
modify this result we will only state the result here of 
how the new value of be affects the published value for 
bn, but this updated value will not be used in the calcula- 
tion of the average value of b np . From the measurement 
of Burgy et al. 72] b c /b H = (-1.756 ± 0.005) fm gives 
b H = (-3.787±0.011) fm (original value: b H = (-3.78± 
0.02) fm. Note that Burgy et al. reported the ratio of 
the atomic scattering lengths as opposed to the nuclear 
scattering lengths: we have converted their ratio into a 
ratio of nuclear scattering lengths using the expression 



(bc/\b H \ato m ) (1 



9ff 



6b r , 



■>H\ 



where b ne = (—1.33 ± 0.03) x 10 -3 fm is the neutron- 
electron scattering length |jjH I M l97 l l9S| . 

To calculate an average of all previous measurements 
we have followed the following procedure given by Hagi- 
wara et al 99 . The data was combined into an aver- 
age with results weighted inversely with the size of their 
(lcr) uncertainties. The result of this evaluation is b np = 
(-3.7410 ± 0.0010) fm and b nd = (6.6727 ± 0.0045) fm. 
The reduced chi square, Xri °f the average value for H is 
3.6 and for D is 0.7. Since the Xr f° r H is greater than 
1 we chose to scale up the uncertainty by yXr based on 
the method of Hagiwara et al [9^ ■ 



We consulted the existing compilations of previous 
measurements [36l l7lj for which we found the measure- 
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